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Abstract 

We uncover geometric aspects that underlie the sum of two independent stochastic 
variables when both are governed by q— Gaussian probability distributions. The 
pertinent discussion is given in terms of random vectors uniformly distributed on a 
p— sphere. 



1 Introduction 



Nonextensive statistical physics provides a rich framework for the interpreta- 
tion of complex systems' behavior whenever classical statistical physics fails 
PQ. The basic tool for this approach is the extension of the classical Boltz- 
mann entropy to the wider class of Tsallis entropies. In this context, the usual 
Gaussian distributions is extended to the q-Gaussian distributions, to be de- 
fined below. The study of the properties of these distributions is an interesting 
problem, being the subject of a number of recent publications pQ. Of special 
interest is the extension of the usual stability result that holds in the Gaussian 
case, namely, that if Xi £ R and X 2 £ M are independent Gaussian random 
variables with unit variance, then the linear combination 



Z = a\X\ + a 2 X 2 



is again Gaussian and 



+ a 2 2 X, 
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where X is Gaussian with unit variance, and ~ denotes equality in distribu- 
tion. 

This stability property is at the core of the central limit theorem (CLT), 
which describes the behavior of systems that result of the additive superpo- 
sition of many independent phenomena. The CLT can be ranked among the 
most important results in probability theory and statistics, and plays an es- 
sential role in several disciplines, notably in statistical mechanics. Pioneers 
like A. de Moivre, P.S. de Laplace, S.D. Poisson, and C.F. Gauss have shown 
that the Gaussian distribution is the attractor of the superposition process of 
independent systems with a finite second moment. Distinguished authors like 
Chebyshev, Markov, Liapounov, Feller, Lindeberg, and Levy have also made 
essential contributions to the CLT-development. As far as physics is concerned 
one can state that, starting from any system with any finite variance distribu- 
tion function (for some measurable quantity x), and combining additively a 
sufficiently large number of such independent systems together, the resultant 
distribution function of x is always Gaussian. 

A natural question is thus the extension of the stability result ([T]) to the nonex- 
tensive case, that is, for q-Gaussian distributions. This interesting problem is 
currently the subject of several publications (see for example [2]) in which 
possible extensions of the CLT to the nonextensive context are studied. The 
aim of this communication is to give some geometric insight into the behavior 
of q-Gaussian distributions for the case q < 1. 



2 Definitions and notations 

In nonextensive statistics, the usual Shannon entropy of a density probability 
fx, namely 

H 1 (X) = - J f x \ogf x 
is replaced by its Tsallis version 

where the nonextensivity index q is a real parameter, usually taken to be pos- 
itive. It can be checked by applying L 'Hospital's rule that Shannon's entropy 
coincides with the limit case 

lim H q (X) = H x (X) 

It is a well-known result that the distribution that maximizes the Shannon 
entropy under a covariance matrix constraint EXX T = K (where K is a 



2 



symmetric definite positive matrix) is the Gaussian distribution 

1 



fx (X) = — -rexp {-XTK^X 

|7TA|2 v 

Its nonextensive counterpart, called a q-Gaussian, is denned as follows. 

Definition 1 The n—variate distribution with zero mean and given covari- 
ance matrix EXX T = K having maximum Tsallis entropy is denoted as 
G q (K) and defined as follows for < q < 1 : 



f x {X) = A il (l-X t E- 1 Xy- 9 , (2) 

with matrix £ = pK , parameter p defined as p = 2jE^ + n an d notation 
(x) + = max (x, 0) . Moreover, the partition function is 

^-r(^)KEr' 

We note that this distribution has bounded support; namely, fx {X) ^ only 
when X belongs the ellipsoid 



£z = {Z eW 1 ; Z ty £T x Z < l} . 



We also need the notion of spherical vector, defined as follows: 

Definition 2 A random vector X G IR n is spherical if its density fx is a 
function of the norm \X\ of X only, namely 

fx(X)=g(\X\) 

for some function g : R + — > R + . 

An alternative characterization of a spherical vector is as follows [3]: 
Proposition 3 A random vector X G M™ is spherical if 

X ~ AX 

for any orthogonal matrix A, where sign ~ denotes equality in distribution. 

This property highlights the importance of spherical vectors in physics since 
they describe systems that are invariant by orthogonal transformation. 

A fundamental property of a spherical vector is the following: 
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Proposition 4 ^ If X G M. n is a spherical random vector, then it has the 
stochastic representation 

X ~rU 

where U is a uniform vector on the sphere S n = jx G M. n ; X T X = 1 j and r is 
a positive scalar random variable independent of U. Moreover, r has stochastic 
representation 

r~\X\. (3) 



3 A heuristic approach 

We start with a heuristic approach to the stability problem, namely the be- 
havior of the random variable Z = a%X\ + a 2 X 2 when X\ and X2 are two unit 
variance, q-Gaussian independent random vectors in M n with nonextensivity 
parameter q < 1; let us assume that the following hypothesis - called (H) 
hypothesis : 

n + —— G N, (4) 
1-q 

holds so that = — 1 where p > n is an integer; a classical result is that 
X\ (resp. X2) can then be considered as the n— dimensional marginal vector 
of a random vector U\ (resp. U 2 ) that is uniformly distributed on the unit 
sphere <S p _i in MP. Thus, there exist random vectors X t and X 2 in W°~ n such 
that 









x 2 






and U2 = 






x K 




x 2 



are two p— dimensional independent vectors uniformly distributed on S p . Then, 
the sum U\ + U 2 is a spherical vector and has stochastic representation 

a\Ui + a 2 U 2 ~ rU 

where U is uniform on S p . Now, by equation ([3]), the random variable r is 
distributed as 

r ~ |aiLq + a 2 U 2 \ = \Ja\ + a 2 + 2\aia 2 
where A = lf[U 2 : this can be easily deduced from 

|aiJ7i + a 2 U 2 \ = yfifiyfUi + a\U%U 2 + 2 ai a 2 UfU 2 

remarking that U\U\ = U 2 U 2 = 1. But A is a random variable with q- 
Gaussian distribution! We prove this result by noticing that, conditioned to 
U 2 = u 2 , random variable A is the angle between JJ\ and the fixed direction u 2 . 
Since U\ is spherical, we may restrict our attention to the angle between Lq 
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and the first vector of the canonical basis in M n , so that we look for the distri- 
bution of the first component of Ui, which follows a q-Gaussian distribution 
with parameter q\ such that 

1 P- 1 i 



1 - qx 



Since this distribution does not depend on our initial choice U 2 — U2, random 
variable A follows unconditionally the above cited distribution. We conclude 
that the n— dimensional marginal Z = a\Xi + a 2 X 2 of vector a-JJ\ + a 2 U 2 is 
distributed as 

a\X\ + a 2 X 2 ~ rX 

where X is the n— dimensional marginal vector of U so that X is again 
q— Gaussian with parameter q. Moreover, this result extends to the case where 
Xi and X 2 both havd"*! a covariance matrix K 7^ / by multiplying vectors Xi 
and X 2 by matrix K2. Consequently, we have deduced the following 

Theorem 5 If Xi and X 2 are two q-Gaussian independent random vectors 
in R n with covariance matrix K and nonextensivity parameter q < 1 and if 
hypothesis (H) holds then 

aiX 1 + a 2 X 2 ~ (ai o a 2 ) X 

where X is again q-Gaussian with same covariance matrix K and same nonex- 
tensive parameter q as X\, and where 



«i «2 = v a? + a 2 + 2Aaia 2 , (5) 

the random variable A being independent of X and again q-Gaussian dis- 
tributed with nonextensive parameter q\ defined by 

= (n- 1) - (n-3)g 
qX (n + 1) - (n- l)q K ) 



Two remarks are of interest at this point: 

• the univariate framework n — 1 is the only case for which random variable 
A has the same nonextensivity parameter q\ as X\ and X 2 ; 

• however, we note that 

lim q\ = 1. 

n— »+oo 

This means that for large dimensional systems, the random variable A con- 
verges to the constant and we recover the deterministic convolution; this 

1 the case where X\ and X 2 have distinct covariance matrices is more difficult and 
left to further study 
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Figure 1. nonextensivity parameter q\ as a function of q for dimensions 
n = 1, 2, 3, 5, 10 and 100 (bottom to top) 

is coherent with the fact that large dimensional q— Gaussian vectors are 
"close" to Gaussian vectors by De-Finetti inequality. 

The curves in Figure [1] show the nonextensive parameter q\ as a function of q 
for several values of dimension n. 

More can be said about the algebra cii o a<i- 



Theorem 6 The algebra a\ o a 2 defined as in (TJP is associative and for any 
n>2, 



n 




As an example, 



ai o a 2 o a 3 = Ja\ + a\ + a| + 2\i 2 aia 2 + 2Ai 3 aia 3 + 2X2^2^. 



PROOF. By definition, 



n 



a 1 oa 2 o...oo B = 





6 



Since \Ui\ = 1, we deduce, by denoting UfUj = Ay, that 



^2 a] + 2^ Xijdidj. 

i=l i<j 

By the same proof as above, we deduce that each A, is q— Gaussian distributed 
with parameter q\. We remark that random variables Ajj are independent 
pairwise but are obviously not mutually independent. 



CL\ o a>2 o 



o a r , 



\ 



4 Generalization 



The preceding result was derived under the hypothesis (H) as expressed by 
(jlj), that is, for specific values of q < 1 only; we show in this section that this 
result holds in fact without this hypothesis - for all values of q < 1 - but the 
proof requires more elaborate analytic tools. Our main result is 

Theorem 7 Theorem^ holds for all values of q such that < q < 1. 



PROOF. The characteristic function associated to the q— Gaussian distribu- 
tion (PJD is 



[\/u T Ku 



2 1 



where Js_i is the Bessel function of the first kind and with parameter | — 1 
where 

p = 2 h n. 

1 - g 

According to Gegenbauer jH 367, eq.16], 

i\ J u (z) Jp (z) = r 

2) Z v z v Jo ' (Z 2 + z 2 -2Zzcos 



rrt v+ 1 -)r( 1 -) J -^l J -^±=[" '^HHH^zl ^ 



Choosing Z = aiVu T Ku, z = a2Vu T Ku, A = — cos0 and v = | — 1, this 
equality can be rewritten as 



where A is distributed according to 



u 



/(A)= j>±i) (i-^r 1 
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Since q\ is denned by 
we deduce (jEJ). 



1 1 _ p- 1 



Let us recall the scaling behavior of Gaussian vectors 



aiXi + a 2 X 2 ~ yaf + c^X 

which can be probabilistically interpreted in the context of a— stable distri- 
butions: a distribution f a is a— stable if, for X\ and X 2 independent with 
distribution f a , the linear combination 

a 1 X 1 + a 2 X 2 ~ (| ai r + |a 2 r)«X, 

where X follows again distribution / a . Thus, a Gaussian distribution is a— stable 
with a = 2. The result of Thm.l can be viewed as follows: q— Gaussians are 
not a— stable (unless q — 1 which corresponds to the Gaussian case a = 2); 
however, their scaling behavior is close to the Gaussian a = 2 case, except for 
the fact that the scaling variable a\ o a 2 includes an additional random term 
A. 



5 Geometric interpretation 



Geometrically, the Gaussian scaling factor W a\ + a 2 can be interpreted, ac- 
cording to Pythagoras' theorem, as the length of the hypotenuse of a right 
triangle with sides of lengths \a\\ and \a 2 \. The q— Gaussian case corresponds 
to a triangle for which the angle between |ax| and | a 2 1 , let us call it (f>, fluctuates 
around rectangularity. 

The distribution of the angle <fi where A = — cos is given by 

f * W = r S t('A 3i °" M ^^'"i^4 



2 \2 



This distributions is shown in Figure [3] for values of the parameter q 
0.99, 0.9, 0.5 and 0.1 (top to bottom). 
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Figure 2. the geometric interpretation of a\ o 02 in the Gaussian case (q = 1 left); in 
the q— Gaussian case (left), a\ o a<i is randomly chosen as one of the hypothenuses 
represented, the angle <f> between sides a% and a2 being distributed as shown on 
Figure [3] 



A 




Figure 3. the distribution of angle <f> for values of the parameter q\ = 0.99, 0.9, 0.5 
and 0.1 (top to bottom). 



We remark that this distribution is symmetric around the angle <f> = ~ and 
that, as q — >• 1, the angle (f> becomes deterministic and equal to |. Further, 
the usual scaling law for Gaussian distributions (CQ) is recovered. 



9 



5.1 An optical analogy 



We remark that formula (jSJ) exhibits a close resemblance with the interference 
formula for the amplitude of the superposition of two optical beams. Interfero- 
metric optical testing is based on these phenomena of interference. Two-beam 
interference is the superposition of two waves, such as the disturbance of the 
surface of a pond by a small rock encountering a similar pattern from a second 
rock. When two wave crests reach the same point simultaneously, the wave 
height is the sum of the two individual waves. Conversely, a wave trough and a 
wave crest reaching a point simultaneously will cancel each other out. Water, 
sound, and light waves all exhibit interference. A light wave can be described 
by its frequency, amplitude, and phase, and the resulting interference pattern 
between two waves depends on these properties, among others. Our present 
interest lies in the two-beam interference equation. It gives the irradiance I 
[S] for monochromatic waves of irradiance 1%, and I 2 in terms of the phase 
difference A lj2 expressed as cos (j) = cos (0i — 2 )- We have 



and, in terms of the A— amplitudes / = A 2 , 

A 2 = A\ + A\ + 1A X A 2 cos (j>. 

If the emission of the two beams could be so arranged that the phase difference 
becomes random [7118119] . this physical analogy would be exact. 

5.2 Study of the composition law o 
The composition law 



has been studied in [5] , in the more general case where a\ and a 2 are indepen- 
dent, positive random variables. The associativity result is as follows 

Theorem 8 [5, p. 18 thm.l] The composition law o is associative if and only 
if either 






or 



or 




(7) 



where X ~ G q (0, 1) for some q > 0. 
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This theorem can be interpreted as follows: the only cases where the compo- 
sition law o is associative is 

(1) the Gaussian case (q — 1, q\ — 1) 

(2) the Cauchy case (g = 2) 

(3) the present q— Gaussian case with < q < 1 

It is thus a remarkable property that the only cases of associativity of this 
composition law correspond to the whole range of q— Gaussian distributions 
with < q < 1 or to the Cauchy case q = 2. We note moreover that in the 
limit case q = 0, a\ o a2 in ([7]) reduces to a Bernoulli random variable that 
takes values a\ + 02 and a\ — 02 with probability \. 

5.3 A Central Limit Theorem for the composition law o 

In the same spirit as the central limit theorem for the usual addition, a central 
limit theorem exists for the composition law o. Before we give its rigourous 
expression as established in [5] , let us look at a special case of it based on the 
theory of superstatistics. Let us consider the random walk 

1 n 

where Xi are independent q— Gaussian random variables with same variance 
a 2 and where q < 1; since all Xi have finite variance, the usual Central Limit 
theorem applies and 

But by theorems §6§ and ([7]), we have also 




with notation 

n 

di = ai o • • • o a n 

i=l 

where X is q— Gaussian with the same nonextensivity parameter q. Since, by 
the superstatistics theory, a Gaussian random variable G can be decomposed 
as 

G = Xp X 

where Xp is chi-distributed with p degrees of freedom, we deduce that the 
following limit 

71 1 " 

i=l a V n a \J n i=i n ^+°° 
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should hold. But this result is easy to check at least under hypothesis (H): in 
this case, 

where 

n 
i=l 

where Ui are independent and uniformly distributed on the sphere <S p _i. By 
the Central Limit Theorem, r n — > \G\ where G is a Gaussian vector in MP; 
hence r n converges indeed to a x distributed random variable with p degrees 
of freedom. 

It turns out that a much more general result holds, namely 

Theorem 9 [5] If {a-i} are positive, independent and identically distributed 
random variables with variance a, then the composition 

1 n 
— i=<g)ai — » x P 

where x P is a chi- distributed random variabl^\with parameter p. 

This result can be considered as a central limit theorem for the algebra o 
defined by (JSj). 



6 Conclusions 



In this work we have uncovered interesting geometric aspects that underlie the 
sum Z of two stochastic variables a±X and 02X2 (ai, 02 are scalars and X±, X2 
are n— variate vectors). The alluded geometry becomes operative when the two 
variables are governed by q— Gaussian probability distributions with q < 1. We 
found that its sum Z turns out to be q— Gaussian with same nonextensivity 
parameter q multiplied by an independent random factor a± o 02. In turn, the 
random factor can be described as a random and symmetric mixture of the two 
constants a\ and 02, the random factor involved following itself a q— Gaussian 
distribution. 



2 we recall that parameter p enters the picture through the distribution of the 
random variable A included in the composition law 
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